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Abstract 

We study continuous time Glauber dynamics for random config- 
urations with local constraints (e.g. proper coloring, Ising and Potts 
models) on finite graphs with n vertices and of bounded degree. We 
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show that the relaxation time (defined as the reciprocal of the spec- 
tral gap |Ai — A2I) for the dynamics on trees and on planar hyper- 
bolic graphs, is polynomial in n. For these hyperbolic graphs, this 
yields a general polynomial sampling algorithm for random configu- 
rations. We then show that for general graphs, if the relaxation time 
T2 satisfies T2 = 0(1), then the correlation coefficient, and the mutual 
information, between any local function (which depends only on the 
configuration in a fixed window) and the boundary conditions, decays 
exponentially in the distance between the window and the boundary. 
For the Ising model on a regular tree, this condition is sharp. 

1 Introduction 

Context 

In recent years, Glauber dynamics on the lattice Z'^ was extensively studied. 
A good account can be found in (25j. In this work, we study this dynamics 
on other graphs. 

The main goal of our work is to determine which geometric properties 
of the underlying graph are most relevant to the mixing rate of the Glauber 
dynamics on particle systems. 

To define a general particle system |.5Tj on an undirected graph G = 
{V, E) , define a configuration as an element a of where A is some finite set, 
and to each edge (f , w) G E, associate a weight function a^w : Ax A IR+. 
The Gibbs distribution assigns every configuration a probability proportional 
to Il{v,w}eE '^v^i'^^^ '^^) ■ The Ising model (for which a^,^„(c^^,, c^^„) = e^'^"'^™) 
and the Potts model are examples of such systems; so is the coloring model 
(for which a^^ = K^^^J 

On a finite graph, the Heat-Bath Glauber dynamics is a continuous time 
Markov chain with the generator 
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where cx^ is the configuration s.t. 
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It is easy to check that this dynamics is reversible with respect to the 
Gibbs measure. An equivalent representation for the Glauber dynamics, 
known as the Graphical representation, is the following: Each vertex has a 
rate 1 Poisson clock attached to it. These Poisson clocks are independent of 
each other. Assume that the clock at v rang at time t and that just before 
time t the configuration was a. Then at time t we replace ct{v) by a random 
spin (j'{v) chosen according to the Gibbs distribution conditional on the rest 
of the configuration: 



We are interested in the rate of convergence of the Glauber dynamics to 
the stationary distribution. Note that this process mixes n = |V^| times faster 
than the corresponding discrete time process, simply because it performs (on 
average) n operations per time unit while the discrete time process performs 
one operation per time unit. 

In section 12.11 we describe a connection between the geometry of a graph 
and the mixing time of Glauber dynamics on it. In particular, we show that 
for balls in hyperbolic tilings, the Glauber dynamics for the Ising model, the 
Potts model and proper coloring with A + 2 colors (where A is the maximal 
degree), have mixing time polynomial in the volume. An example of such a 
hyperbolic graph can be obtained from the binary tree by adding horizontal 
edges across levels; another example is given in Figure ^ 

In sections I2.3I E1 we study Glauber dynamics for the Ising model on reg- 
ular trees. For these trees we show that the mixing time is polynomial at 
all temperatures, and we characterize the range of temperatures for which 
the spectral gap is bounded away from zero. Thus, the notion that the two 
sides of the phase transition (high versus low temperatures) should corre- 
spond to polynomial versus super-polynomial mixing times for the associ- 
ated dynamics, fails for the Ising model on trees: here the two sides of the 
high/intermediate versus low temperature phase transition just correspond 
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Figure 1: A ball in hyperbolic tiling 



to uniformly bounded versus unbounded inverse spectral gap. We also exhibit 
another surprising phenomenon: On infinite regular trees, there is a range 
of temperatures in which the inverse spectral gap is bounded, even though 
there are many different Gibbs measures. 

In section El of the paper we go beyond trees and hyperbolic graphs and 
study Glauber dynamics for families of finite graphs of bounded degree. We 
show that if the inverse spectral gap of the Glauber dynamics on the ball 
centered at p stays bounded as the ball grows, then the correlation between 
the state of a vertex p and the states of vertices at distance r from p, must 
decay exponentially in r. 



The graphs. Let G = (V, E) be an infinite graph with maximal degree A. 
Let p be a distinguished vertex and denote by Gr = (K, Ej.) the induced 
graph on Vr = {v E V : dist(p, f) < r}. Let Ur be the number of vertices in 
Gr- At some parts of the paper we will focus on the case where G = T = 
{V, E) is the infinite 6-ary tree. In these cases, Tr^'' = Er) will denote 
the r-level 6-ary tree. 

The Ising model. In the Ising model on Gr at inverse temperature (3, every 
configuration a G {—1,1}^'' is assigned probability 
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where Z{P) is a normalizing constant. When Gr = Tr , this measure has 
the following equivalent definition [H]: Fix e = (1 + e^^)~^. Pick a random 
spin ±1 uniformly for the root of the tree. Scan the tree top-down, assigning 
vertex v a spin equal to the spin of its parent with probability 1 — e and 
opposite with probability e. 

The Heat-Bath Glauber dynamics for the Ising model chooses the new 
spin (j'{v) in such a way that: 

PW(v) = +l\a] / ^ , x\ 

^ ^ = exp 2/5 > a H • 

P a'(v) = -1 a ^ V ^ ^ ^ V 

See [21] or j2Sl for more background. 



Mixing times. 

Definition 1.1. For a reversible continuous time Markov chain, let = 
^1 < A2 < . . . < Afc be the eigenvalues of —C where C is the generator. The 
spectral gap of the chain is defined as X2, and the relaxation time, T2, is 

defined as the inverse of the spectral gap. 

Note that the corresponding discrete time Glauber dynamics has transi- 
tion matrix M = I + -C, where n is the number of vertices. Moreover, the 
eigenvalues of M are 1,1 — — ,1 — and therefore the spectral gap of 

the discrete dynamics is the spectral gap of the continuous dynamics divided 
by n. 

Definition 1.2. For measures /i and v on the same discrete space, the total- 
variation distance, dv{lJ',i^), between fi and v is defined as 

X 

Definition 1.3. Consider an ergodic Markov chain {Xt} with stationary 
distribution n on a finite state space. Denote by the law of Xt given 
Xq = X. The mixing time of the chain, ti, is defined as 

n = inf{t : supdy(P*,P*) < e"^}. 

x,y 
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For t > ln(l/e)ri, we have 

suprfy(P^,7r) < suprfy(P*.,P*) < e. 

X x,y 

Using T2 one can bound the mixing time ri, since every reversible Markov 
chain with stationary distribution vr satisfies (see, e.g., P), 

T2<n< T2 (l + log ((min7r(a))-i)) . (2) 

For the Markov chains studied in this paper, this gives T2 < ti < 0{n)T2. 

Cut- Width and relaxation time. 

Definition 1.4. The cut-width ^(G) of a graph G is the smallest integer 
such that there exists a labeling fi, . . . ,f„ of the vertices such that for all 
1 < k < n, the number of edges from {vi, . . . ,Vk} to {ffc+i, • • • ,Vn}, is at 
most C,{G). 

Remark: The vertex-separation of a graph G is defined analogously to 
the cut-width in terms of vertices among {vi, . . . ,Vk} that are adjacent to 
{vk+i, . . . ,Vn}- In j2ni it is shown that the vertex-separation of G equals its 
path-width, see 1221 ■ [13 the cut- width was called the exposure. 

Generalizing an argument in fl^ Theorem 6.4] for Z'^, (see also we 
prove: 

Proposition 1.1. Let G be a finite graph with n vertices and maximal degree 
A. 

1. Consider the Ising model on G. The relaxation time of the Glauber 
dynamics is at most ne^^^^'^^'^'^^^^ . 

2. Consider the coloring model on G. If the number of colors q satisfies 
g > A + 2, then the relaxation time of the Glauber dynamics is at most 
(A + l)n(g-l)«(^)+i. 

Analogous results hold for the independent set and hard core models. 
Cut- Width and long-range correlations for hyperbolic graphs. The 

usefulness of Proposition ^3 comes about when we bound the relaxation time 
of certain graphs by estimating their cut-width. The following proposition 
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bounds the cut-width of balls in hyperbolic tilings of the plane. Recall that 
the Cheeger constant of an infinite graph G is 

A C G;0 < |A| < ooj , (3) 

where dA is the set of vertices of A which have neighbors in G\A. 

Proposition 1.2. For every c > and A < oo, there exists a constant 
C = C{c, A) such that if G is an infinite planar graph with 

• Cheeger constant at least c, 

• maximum degree bounded by A and 

• for every r no cycle from Gr separates two vertices of G\ Gr, 

then C,{Gr) < Clogn,, for all r, where is the number of vertices of Gr- 

Combining this with Proposition 11.11 we get that the Glauber dynamics 
for the Ising models on balls in the hyperbolic tiling has relaxation time 
polynomial in the volume for every temperature. On the other hand, we 
have the following proposition: 

Proposition 1.3. Let G be a planar graph with bounded degrees, bounded co- 
degrees and a positive Cheeger constant. Then there exist [3' < oo and 5 > 
such that for all r, all j3 > j3' , and all vertices u,v in Gr, the Ising model 
on Gr satisfies that E[(Ttj(T^] > 6. In other words, at low enough temperature 
there are long-range correlations. 

This shows that for the Ising model on balls of hyperbolic tilings at very 
low temperature, there are long-range correlations coexisting with polyno- 
mial time mixing. While there is no characterization of all planar graphs 
with positive Cheeger constants, an important family of such graphs is Cay- 
ley graphs of nonelementary fuchsian groups which are nonamenable, and 
hence have a positive Cheeger constant, see [SH 1211 UH] for background and 
[TT] for some explicit estimates. 

Relaxation time for the Ising model on the tree. The Ising model 
on the 6-ary tree has three different regimes, see E]- In the high tem- 
perature regime, where 1 — 2e < 1/6, there is a unique Gibbs measure on 
the infinite tree, and the expected value of the spin at the root dp given any 
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boundary conditions (JgrpW decays exponentially in r. In the intermediate 

regime, where 1/6<1 — 2e< 1/Vb, the exponential decay described above 
still holds for typical boundary conditions, but not for certain exceptional 
boundary conditions, such as the all + boundary; consequently, there are 
infinitely many Gibbs measures on the infinite tree. In the low temperature 
regime, where 1 — 2e > l/Vb, typical boundary conditions impose bias on 
the expected value of the spin at the root cxp. 

Theorem 1.4. Consider the Ising model on the b-ary tree Tj. = Tr^^ with r 
levels. Let e = (1 + e^^)'"^ . The relaxation time T2 for Glauber dynamics on 
can be bounded as follows: 

1. The relaxation time is polynomial at all temperatures: T2 = n?^^"^^^^^^\ 
Furthermore, the limit 

r^oo log(nr) 

exists. 

2. Low temperature regime: 

(a) If 1 — 2e > 1/Vb then snprT2{Tr) = 00. In fact, T2{Tr) = 
nirir^"^''^'^'^'^'^) when 1 - 2e > 1/Vb and r2(T,) = l](logn,) when 
1 - 2e = 1/y/b. 

(b) Moreover, the degree of T2 tends to infinity as e tends to zero: 

,rp\ n(log(l/e)) 

3. Intermediate and high temperature regimes: 

// 1 — 2e < 1/ Vb then the relaxation time is uniformly bounded: T2 = 
0(1) . Furthermore, this result holds for every external field {H{v)}y^Tr ■ 

In particular we obtain from Equation Q that in the low temperature 
region ti = nf^^\ and in the intermediate and high temperature regions 
Ti = 0{nr). 

A recent work by Peres and Winkler |i33j compares the mixing times of 
single site and block dynamics for the heat-bath Glauber dynamics for the 
Ising model. 

They show that if the blocks are of bounded volume, then the same mixing 
time up to constants is obtained for the single site and block dynamics. 
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high 


< 1/2 


unbiased 





0(1) 


med. 




biased 


^ 


0(1) 


low 




biased 


inf> 




freeze 


l-o(l) 


biased 


l-o(l) 





Table 1: The Ising model on binary trees. Here the root is denoted p, and 
the vertices at distance r from the root are denoted dT. 

Combining these results with the path coupling argument of Section HI it 
follows that Ti = 0(logn,.) in the intermediate and high temperature regions. 

We emphasize that Theorem 11.41 implies that in the intermediate region 
1/2 < 1 — 2e < 1/v^, the relaxation time is bounded by a constant, yet, in 
the infinite volume there are infinitely many Gibbs measures. This Theorem 
is perhaps easiest to appreciate when compared to other results on the Gibbs 
distribution for the Ising model on binary trees, summarized in Table Q 

The proof of the low temperature result is quite general and applies to 
other models with "soft" constraints, such as Potts models on the tree . 
Spectral gap and correlations. At infinite temperature, where distinct 
vertices are independent, the Glauber dynamics on a graph of n vertices 
reduces to an (accelerated by a factor of n) random walk on a discrete n- 
dimensional cube, where it is well known that the relaxation time is 6(1). 
Our next result shows that at any temperature where such fast relaxation 
takes place, a strong form of independence holds. This is well known in Z"^, 
see (23, but our formulation is valid for any graph of bounded degree. 

Theorem 1.5. Denote by the configuration on all vertices at distance 
r from p. If G has bounded degree and the relaxation time of the Glauber 
dynamics satisfies T2{Gr) = 0(1), then the Gibbs distribution on Gr has the 
following property. For any fixed finite set of vertices A, there exists ca > 
such that for r large enough 

Cov[f,g] < e-'^- VVar(/)Var(^7) , (4) 

provided that f{a) depends only on a a and g{cr) depends only on o^. Equiv- 
alently, there exists dj^ > such that 

I[aA,cTr]<e-''^\ (5) 
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where I denotes mutual information (see J^.) 

This theorem holds in a very general setting which includes Potts models, 
random colorings, and other local-interaction models. 

Our proof of Theorem ll.5l uses "disagreement percolation" and a coupling 
argument exploited by van den Berg, see j2] , to establish uniqueness of Gibbs 
measures in Z"^; according to F. Martinelli (personal communication) this 
kind of argument is originally due to B. Zegarlinski. Note however, that 
Theorem 11.51 holds also when there are multiple Gibbs measures - as the 
case of the Ising model in the intermediate regime demonstrates. Moreover, 
combining Theorem 11.51 and Theorem II. 41 one infers that for 1 — 2e < 1/Vb, 
we have limr^oo -^[^o, cTt] = 0. This yields another proof of this fact which 
was proven before in jHl [T21 El • 

Plan of the paper 

In section|2]we prove Proposition ll. ll via a canonical path argument, and give 
the resulting polynomial time upper bound of Theorem 11.41 part ^ We also 
present a more elementary proof of the upper bound on the relaxation time 
for the tree, which gives sharper exponents and the existence of a limiting 
exponent; this proof uses Martinelli's block dynamics to show sub-additivity. 
In section |3] we sketch a proof of Theorem 11.41 part and present a proof 
of Theorem 11.41 part I2bl These lower bounds are obtained by finding a 
low conductance "cut" of the configuration space, using global majority of 
the boundary spins for the former result, and recursive majority for the 
latter result. In section |3 we establish the high temperature result, using 
comparison to block dynamics which are analyzed via path-coupling. Finally, 
in section El we prove Theorem 11.51 bv a Peierls argument controlling "paths 
of disagreement" between two coupled dynamics. 

Remark: Most of the results proved here were presented (along with 
proof sketches) in the extended abstract [T^. However, the proofs of our 
results for hyperbolic graphs (see Section 2.2), which involve some interesting 
geometry, were not even sketched there. Also, the general polynomial upper 
bound for trees that we establish in Section 2.3 is a substantial improvement 
on the results of jTHj, since it only assumes the dynamics is ergodic and allows 
for arbitrary hard-core constraints. 
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2 Polynomial Upper Bounds 



2.1 Cut- Width and mixing time 

We begin by showing how part ^ of Proposition II . II imphes the upper bound 
in part ^ of Theorem 11.41 

Lemma 2.1. Let tJ^^ be the b-ary tree with r levels. Then, ^(tJ'''*) < (6 — 
l)r + 1. 

Proof. Order the vertices using the Depth first search left: to right order , i.e., 
use the following labeling for the vertices: The root is labeled (0, 0, . . . , 0). 
The children of the root are labeled (1,0,..., 0) through (6, 0, . . . , 0), and so 
on, so that the children of (oi, 02, . . . , a^, 0, . . . , 0) are (oi, 02, . . . , flfc, 1, . . . , 0) 
through (ai, 02, . . . , a^, 6, . . . , 0). Then order the vertices lexicographically. 
Note that in the lexicographic ordering, a vertex always appears before its 
children. When we enumerated all vertices up to (ai, 02, . . . , flr), the only 
vertices that were enumerated but whose children were not enumerated are 
among the set of at most r vertices 

{(0, 0, . . . , 0), (ai, 0, . . . , 0), (ai, 02, . . . , 0), . . . , (oi, 02, ... , a^)} . 

Each of these vertices has at most b children, and for all but (ai, 02, ... , ctr) 
at least one child has already been enumerated. Therefore, 

e(r(^))<(&-i)r + i. 

□ 

Corollary 2.2. 1. The relaxation time of the Glauber dynamics for the 
Ising model on Tr''^ is at most 

2. The relaxation time of the Glauber dynamics for the coloring on Tr^^ 
with q > b + 2 colors is at most 

(6 + l)nJ+2(''-i)i°s<>(«) 
Proof. The Corollary follows from Lemma 12.11 and Proposition 11.11 □ 
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The upper bound in part of Theorem 11.41 follows immediately. 



Proof of partUl of Proposition \l.A The proof follows the lines of the proof 
given in [23 Theorem 6.4] for the Ising model in Z'^, (see also [T3]). 

Let r be the graph corresponding to the transitions of the Glauber dy- 
namics on the graph G. Between any two configurations a and t], we define 
a "canonical path" 7(0",?]) as follows. Fix an order < on the vertices of G 
which achieves the cut-width. Consider the vertices < f 2 < . . . at which 
ay ^ Tjy. 

We define the fc-th configuration a'''^^ on the path 'j(a,r]) by giving spin 
(Ty to every labeled vertex v < v^, spin rj^ to every labeled vertex v > Vk, 
and spin ay = rjy for every unlabeled vertex v. Note that a^^^ = t] and 
0-id(<^^v)) = (j_ Since a^''~^^ and o"'-^^ are identical except for the spin of vertex 
Vk, they are adjacent in T. This defines 'y{a,ri) (see FigureE)). Note that for 
every k, there are at most C,{G) pairs of adjacent vertices {vi,Vj) such that 
i < k < j, hence any configuration on the canonical path between a and rj 
will have at most ^{G) edges between spins copied from a and spins copied 
from 1]. 

Using canonical paths to bound the mixing rate. For each directed 
edge e = {uj,C) on the configuration graph F, we say that e G '~f{a,ri) if lj 
and C are adjacent configurations in 7(0", 77) . Let 
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Figure 2: The canonical path from a to rj. The vertices on which a and rj 
agree are marked in grey; the other vertices are colored black if their spin is 
chosen according to a and white if their spin is chosen according to rj. 



12 



00000100 
00 »0#00«^ 



• 0«#0000 



•OtlOlOO*" 

00 00000 On 



Figure 3: The injection from (e, to (cr, r/). The vertices on which both 
endpoints of e and ip agree are marked in grey; the other vertices are colored 
black if they precede f^o ^"^^ their spin is chosen according to ip, or if they 
are preceded by and their spin is chosen according to the endpoints of e; 
and are colored white otherwise. 

where /i is the stationary measure (i.e. the Gibbs distribution), and for any 
two adjacent configurations u and C, Q(e) = Q((a;, C)) = //[cuji^fti; — > Q. If L 
is the maximal length of a canonical path, then by the argument in 15;, '2")| , 
the relaxation time of the Markov chain is at most 



Since L < n, it follows that T2 < np, thus it only remains to prove an upper 
bound on p. 

Analysis of the canonical path. For each directed edge e in F, we 
define an injection from canonical paths going through e in the specified 
direction, to configurations on G. To a canonical path 'y{a,T]) going through 
e, such that e = (a'^^^^-', cr^^-'), we associate the configuration ip which has 
spin rjy. for every Vi s.t. i > k and spin ay. for every Vi s.t. i < k. To 
verify that this is an injection, note that one can reconstruct a and rj by first 
identifying the unique ko s.t. iv and ( differ on Vk^ and then taking (as in 
Figure EI) 



T2 < Lp. 



(6) 




and 
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By the property of our labeling, 

^M^M < ^^[a^'-'^]^^W^^^^^ . (7) 

and K > exp(— 2A/3). Now a short calculation concludes the 

proof: 



p < sup ^ 



(7,»? S.t. ee7{(7,»y) 



The last inequality follows from the fact that the map 7 — >■ 95 is injective and 
therefore < 1- □ 

Proof of part of Proposition \l.ll The previous argument does not directly 
extend to coloring, because the configurations cr*-*^-' along the path (as defined 
above) may not be proper colorings. Assume that g > A + 2 and let vi < 
^2 ■ ■ ■ < be an ordering of the vertices of G which achieves the cut-width. 
We construct a path 7(0", rj) such that 

|7(a,77)| < (A + IK (10) 

Moreover, for all r G 7(0", 77) there exists a k such that 

" \ if V > Vk and f 7^{f , . . . , Wfc} 

The way to construct a path 7(0", rj) satisfying (fTTH) and (fTT|l is the following: 
(T° = a. Given a'', we proceed to create a^^^ as follows: Let i{k) = inf{j : 

(T^ if 7^ fi(fc) 
r]^ if = fi(fc) 

is a legal configuration, then a''^^ = p. otherwise, let 

h{k) = mi{j : a^^. = r]^^^^^ and ~ fi(fc)}, 

14 



and let c be a color that is different from ?7^.(j.j and is legal for Vh(k) under a^. 
Such a color exists because g > A + 2. Then, we take 

\ c \iv = Vh(k) 

It is easy to verify that the path satisfies (fTUI) and (fTT|) . Since all legal 
configurations have the same weight, ((Tj) is replaced by 

^^[a]^,U = f^[a^'-'^M^] (12) 

On the other hand, the map 7 ^ (/? is not injective. Instead, by (jll|) . there 
are at most {q — lY^^^ paths which are mapped to the same coloring. We 
therefore obtain that for the coloring model p < n(q — l)?(*^)+i and therefore 
from (Uni) and (ini), 

r2 < (A + l)nq{q - l)^^^). 

□ 

2.2 Hyperbolic graphs 

In this subsection we show that balls in a hyperbolic tiling have logarithmic 
cut-width. Let G = {V, E) be an infinite planar graph and let o & V . Let Gr 
be the ball of radius r in G around o, with the induced edges. The following 
proposition implies Propositions 11.21 and 11.31 

Proposition 2.3. 1. Suppose G has 

• a positive Cheeger constant, 

• degrees bounded by A, 

• for all r, no cycle from Gr separates two vertices of G\ Gr, 

then there exists constants ai and a2 s.t. C,{Gr) < aiA logdC,.!) +q;2A. 

2. Assume that G has bounded degrees, bounded co-degrees, no cycle from 
Gr separates two vertices ofG\Gr, and the following weak isoperimetric 
condition holds: 

\dA\>G\ogi\A\) (13) 
for every finite A (1 G and for some constant C. 
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Then there exist (3' < oo and 5 > Q s.t. for every j3 > j3' , for every r 
and for every u, v in Gr, the free Gibbs measure for the Ising model on 
Gr with inverse temperature [3 satisfies cov(cr„, o"^,) > 5. 

Proof of partU\ of Proposition \2. Consider a planar embedding of G. Since 
no cycle from G^ separates two vertices of G \ G^; all vertices of G \ G,. are 
in the same face of G^, and without loss of generality we can assume that it 
is the infinite face of our chosen embedding of Gr- 

Let T be a shortest path tree from o in G^. In other words, T is a tree 
such that for every vertex f , the path from o to f in T is a shortest path in 
Gr- Let Ci G T be an edge adjacent to o. We perform a depth-first- search 
traversal of T, starting from a = vq, traversing ci to its end vertex vi, and 
continuing in counterclockwise order around T. This defines a linear ordering 
vo ^ vi < ■ ■ ■ < Vn-i of the vertices of Gr- 

Consider the induced ordering wi < W2 < ■ ■ ■ < Wk on the vertices of Gr 
which are at distance exactly r from o. 

Fix i < j- We first consider edges between 

Vij = {u E Gr : Wi < u < Wj,u not ancestor of Wj in T} 

and Gr \ Vij- Note that Vij does not contain any vertex on the paths in T 
from to either Wi or Wj. Obviously, there can be edges from Vij to vertices 
on the paths from a to Wi or Wj. Let e = {u, f } be an edge with one endpoint 
in Vij and the other end in G \ G^ but not on Path(wj) or Path(wj), where 
Path(wj) denote the path in T from o to Wj- Without loss of generality, 
assume that Wi < u < Wj < v- The case where v < Wi < u < Wj is treated 
similarly. 

The path from o to u in T, followed by the edge e, followed by the path 
from t> to o in T, defines a cycle Ge in Gr (see FigureEj). Since Wi < u < Wj < 
V, Ge must enclose exactly one of Wj and Wi- Since the graph is embedded 
in the plane, the ones among those cycles which enclose Wj must form a 
nested sequence, and therefore there is an outermost such cycle Ge* with 
an associated "outermost" edge e* = {u*,v*}. Similarly, among the edges 
such that the corresponding cycle encloses Wi, there is an "outermost" edge 
f* = {x\y*}- 

There can only be edges from Vij to the vertices enclosed by Ge* or by 
Gf* (note that this includes the paths from o to Wi and to Wj) . Since all the 
vertices of G \ G^ are in the infinite face of Gr, hence outside Ge*, the set 
of vertices enclosed by Ge* is the same in G as in Gr- Let A denote the set 
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Figure 4: The cycle Cg defined by e. 

of vertices enclosed by Ce* (including Ce*). We have: \dA\ < 2r + 1, hence 
1^1 < (2r + l)/c, where c is the Cheeger constant of G. Reasoning similarly 
for Cf*, we obtain that the set of vertices in Gr \ Vij adjacent to Vij has size 
at most (4r + 2) /c. 

Let Bj = Vj-ij for j > 2, and Bi = {u ^ Gr '■ u < Wi or u > 
Wk,u not an ancestor of Wi in T}. Let us bound the cardinality of Bj. As 
above, we define Cg* and Cf*. Let A denote the union of Bj, of the vertices 
enclosed by Cg*, and of the vertices enclosed by Cf*. 

Since the vertices of Bj arc at distance at most r — 1 from o, they have 
no neighbors in G \ Gr- Thus the neighborhood of A in G is such that 
dA C Ce- U Cf*, hence \dA\ < 4r + 2, and so \Bj\ < \A\ < (4r + 2)/c. 

Finally, to compute the cut- width, let S = {u : Vq < u < Vi}, and let j 
be such that Wj-i < Vi < Wj. We have: 

Vij-i QS CBiU Vij_i U Bj U Path(wi) U Path(wj_i) U Path(wj). 

Thus the set of edges between S and Gr\S has size at most A(4r + 2)/c + 
(1^1 U Bj\)A + (3r + 1)A, which is at most (3r + 1 + (12r + 6)/c)A. 
Since G has positive Cheeger constant, 

= \Gr-l U dGr-l\ > \Gr-l\{l + c), 

and so \Gr\ > (1 + c)'', that is, r < log |Gr|/ log(c+ 1). Hence the set of edges 
between 5" and G\ 5" has size at most (3 + 12/c)(A/ log(c+ 1)) log |G^| + (1 + 
6/c) A. This concludes the proof. □ 
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Figure 5: The region between Path(M/j) and Path(iyj) 

Proof of partlE of Proposition \2.3i . We use the Random Cluster representa- 
tion of the Ising model (see, e.g. and a standard Peierls path-counting 
argument. For every u and v in Gr, cov((T„, a,,) is the probability that u is 
connected to v in the Random Cluster model. Fix p < 1. The exact value of 
p will be specified later. 

Then, /3 is large enough, i.e., if (1 — e~^)/e~^ > 2p/(l — p), then the 
Random Cluster model dominates percolation with parameter p. So, what 
we need to show is that for a graph satisfying the requirements of part El of 
the proposition and p high enough, there exists 6 > s.t. for every r and 
every u, v in Gr, we have Pp{u ^ v) > 6. By the FKG inequality (see pUj). 

Pp{u ^ v) > Pp{u ^ o)Pp{v o) 

where o is the center. Therefore we need to show that P{v ^ o) is bounded 
away from zero. To this end, we will pursue a standard path counting tech- 
nique: in order for o and v not to be connected, there needs to be a closed 
path in the dual graph that separates o and v. 

Claim 2.4. There exists M = M{G) s.t. for every r and v & Gr there are 
at most M'^ paths of length k in the dual graph of Gr that separate a from v. 

By Claim|231 if we take p > 1 — 1/ (2M) and choose f3 accordingly then the 
probability that there exists a closed path in the dual graph that separates 
and V is bounded away from 1. □ 



Proof of Claim Here again, we consider an embedding of Gr such that 



all the vertices oi G\Gr lie on the infinite face F of Gr ■ 

Let 7 be a shortest path connecting t> to o in G^- Every dual path 
separating v from a must intersect 7. For an edge e let Afc(e) be the set of 
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dual paths of length k separating o from v such that ip intersects e. If A 
is the maximal co-degree in G then |Afc(e)| < A'^ for every e. 

Let e G 7 be such that d{e, o) > exp{k/C) + kA, d{e, v) > exp{k/C) + kA, 
and d{e,F) > kA. We will now show that |Afc(e)| = 0. Assume, for a 
contradiction, that ip G Afc(e). Since ip has length k and d{e,F) > Ak, ip 
does not touch the outer face F, and so the area enclosed by ip in Gr equals 
the area enclosed by ip in G. The dual path ip encloses either v or o. Without 
loss of generality, assume that it encloses o. Let e' be the edge of 7 closest to 
o which if) intersects. Since ip has length k, we get that d{e\o) > exp{k/C), 
and so at least 1 + exp(fc/C) vertices of 7 are enclosed by ip. By (fT^ this 
implies that ip has length strictly greater than k, a contradiction. 

Thus, the total number of paths of length k separating o from v is at 
most 

J2 \Me) I < [2 exp(A;/C) + kA) + AA;] A^ 

e:|Afe(e)|7tO 

□ 

Remark: An isoperimetric inequality of the type of (jl3|) is necessary. An 
example where all other conditions of part El of Proposition 12.31 are satisfied 
and yet the conclusion does not hold can be found in Figure IHl 




Figure 6: An example of a planar graph s.t. for all temperatures, correlations 
decay exponentially with distance. 

2.3 A polynomial upper bound for trees 

In this subsection we give an improved bound on relaxation time for the tree. 
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Let A be a finite set, and let a^^ : A ^ R+ be a weight function. 
Let G be a grapli. Let tfie Glauber dynamics be as defined above, and let 
L = C{A, a, G) be its generator. We say that the Glauber dynamics on 
{A,a,G) is ergodic if for every two legal configurations ai and (T2, we have 
(exp(£))^^^^ > 0. We will prove the following proposition: 

Proposition 2.5. Let b > 2, and let T denote the infinite b-ary tree, and let 
Tn be the b-ary tree with n levels. If the Glauber dynamics on {A, a, T„) is 
ergodic for every n then 

limsup - log {t2 {C {A, a,T„))) < 00. 

Conjecture 2.6. Let b > 2, T denote the infinite b-ary tree, and let Tn be 
the b-ary tree with n levels. If the Glauber dynamics on (A, a,T) is ergodic 
then there exists < r < 00 s.t. 

hm - log (r2(/: (A, a, T„)))=r. (14) 

n— >oo n 

We prove a special case of Conjecture 12.61 

Proposition 2.7. If the interactions are soft, i.e. a^^(a, 6) > for all v, w, a 
and b, then [T^ holds. 

The main tool we use for proving Propositions 12.51 and 12.71 is block dy- 
namics (see e.g. For a spin (or a color) a E A, we denote by £(a, a, n) 
the Glauber dynamics on the 6-ary tree of depth n, under the interaction ma- 
trix a and with the boundary condition that the root has a parent colored 
a. With a slight abuse of notations, we say that T2{a,a,n) is the relaxation 
time for £(a, a, n). 

Lemma 2.8. Let 

f2{a, n) = sup T2{a, a, n) 

Then, for all m and n, 

f2(a, n + m) < hia, n)f2{a, m). 

Proof. Let I = n + m. Partition the tree T/ into disjoint sets Vi,..., 14 to 
be specified below. We call Vi, blocks, and consider the following block 
dynamics: Each block Vi has a (rate 1) Poisson clock, and whenever it rings, 
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Vi updates according to its Gibbs measure determined by the boundary con- 
ditions given by the configurations of t/^'' — Vi and by the external boundary 
conditions. We denote by = C^{Vi, V^) the generator for the block 
dynamics, and let £f be the generator for the block dynamics with the 
boundary condition that the parent of the root has color a. 
By [221 [Proposition 3.4, page 119], 

hia, I) < sup f2(a, Vi) ■ sup r2(£f ) 

i a^A 

We now define the partition to blocks. For every vertex v up to depth n, the 
singleton {f } is a block, and for every vertex w at depth n, the full subtree 
of depth m starting at w is a block (see Figure [7j). All we need now to finish 




Figure 7: Partition of a tree to blocks 



the proof is the following easy claim: 
Claim 2.9. 

supr2(£f ) = f2{a,n). 

Proof. We use the following fact (that could also serve as a definition of 
the relaxation time). Given the dynamics C we define the Dirichlet form 
£[9,9] = lJ:.,rf^m[a - T]igia)-gir))\ Then 



T2 = sup 



£[9,9] 



: fi[g] = . 



(15) 



Clearly, the expression in (fT3j) evaluated for / and £f is equal to the one 
evaluated for g and C{a, a, n), if 



g{ri) = /(cr) for all t] and a s.t. t] \t„ = <J. 



(16) 



Therefore, we need to show that the maximum in (fT3j) for the dynamics £f is 
obtained at a function that satisfies (jl6|) . The maximum in (jl5p is obtained 
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at an eigenfunction of £f . Moreover for every function g, C^{g) satisfies 
with some function /. It now follows that the maximum is obtained at 
a function that satisfies (fTBj) . □ 

□ 

Proof of Proposition \2.^ From Lemma 12.81 and the sub-additivity lemma, 
we learn that ^ 

limsup - log (f2 {C {A, a, T„))) < oo 

By another application of Matinelli's block dynamics lemma, we get that 

T2 {C (A «, T„)) < r2 {C (A a, T,)) ■ {C {A, «, T„_i)) (17) 

and the proposition follows. □ 

Proof of proposition "WTi From Lemma ITHl and the sub-additivity lemma, we 
learn that there exists < r < cx) s.t. 

lim - log {f2 {C {A, a, T„))) = r. 

n^oo n 

For every a, let /ia be the Gibbs measure for the tree of depth n with the 
boundary condition that the parent of the root has color a. Note that /ia is 
the stationary distribution of £(a, a, n). Since the interactions are soft, there 
exists < C < oo s.t. for every a, every n, and every two configurations a 
and rj on the tree of depth n, 

and ^ 

—C{a,n)„^ri < C{a,a,n)„^r, < CC{a,n)„^r,- 

Therefore, by (fT^. 

i^r2 (£ (A a, r„)) < f2 {c (A «, r„)) < (£ (A «, 

and the proposition follows. □ 
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3 Lower Bounds 

Proof of Theorem \1.4[ part\2al Theorem 11.41 part I2alis a direct consequence 
of the extremal characterization of T2 given in ()15p. apphed to the particular 
test function g which sums the spins on the boundary of the tree. It is easy 
to see that fi[g] = and that 



S[g,g]<J2f^[^]KW^^]=0{ 
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We repeat the variance calculation from [8]. When 6(1 — 2e)^ > 1: 



= If ■ 1^1 + ^(6- 1)6^-^(1 -2e) 
It now follows by ^ that if 6(1 - 2e)2 > 1 then 

T, = n (^nlog.(Ml-2e)2)^ 

as needed. Repeating the calculation for the case 6(1 — 2e)^ = 1 yields that 

T2 = fi(logn^). 

The proof follows. □ 

Remark: Suppose that /i admits a Markovian representation where the 
conditional distribution of ctu given its parent o"^, is given by an |^| x |^| 
mutation matrix P. Let A2(-P) be the second eigen- value of P (in absolute 
value), and x the corresponding eigen- vector, so that Px* = A2(-P)x* and 

\x\2 = 1. 

Let g be the test function g = CnX^, where c„(z) is the number of boundary 
nodes that are labeled by i. It is then easy to see once again that S[g,g] = 
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0{nj.). Repeating the calculation from it follows that if 6|A2(-P)P > 1, 
then 



In order to prove the lower bound on the relaxation time for very low tem- 
peratures stated in Theorem II .41 part l2b | we apply (fT5|l to the test function g 
which is obtained by applying recursive majority to the boundary spins; see 
[2HI for background regarding the recursive- majority function for the Ising 
model on the tree. For simplicity we consider first the ternary tree T, see 
Figure 5. Recursive majority is defined on the configuration space as follows. 
Given a configuration a, first label each boundary vertex v by its spin a^- 
Next, inductively label each interior vertex w with the label of the majority 
of the children of w. The value of the recursive majority function g is then 
the label of the root. We write for the spin at v and m„ for the recursive 
majority value at v. 

Lemma 3.1. If u and w are children of the same parent v, then P[m„ 7^ 
< 2e + Se^. 

Proof: 

P[m„ 7^ m^] < P[cr„ 7^ m„] + P[cr^ 7^ m^] + P[(T„ 7^ cr^] + P[cr^ 7^ a^]. 

We will show that recursive majority is highly correlated with spin, i.e. if e 
is small enough (say e < 0.01), then P[m^, 7^ a^] < 4e^. 

The proof is by induction on the distance i from v to the boundary of 
the tree. For a vertex v at distance i from the boundary of the tree, write 
Pi = P[m^ 7^ cTy]. By definition Po = < 4e^. 




Thus in this case 





Figure 8: The recursive majority function. 
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For the induction step, note that if 7^ then one of the following 
events hold: 

• At least 2 of the children of have different cr value than that of a^, 



• One of the children of v has a spin different from the spin at w , and for 
some other child w we have 7^ cr^, or 

• For at least 2 of the children of f , we have 7^ frtw- 

Summing up the probabilities of these events, we see that pi < 3e^ + Gepi^i + 
3pj_i^. It follows that pe < 4e^, hence the Lemma. □ 

Proof of Theorem \1.4\ part IsR Let m be the recursive majority function. Then 
by symmetry E[m] = 0, and E[m^] = 1. By plugging m in definition ()15|). 
we see that 



Observe that if a, r are adjacent configurations (i.e., P[cr — r] > 0) such 
that m(cr) = 1 and m(r) = —1, then there is a unique vertex Vr on the 
boundary of the tree where a and r differ. Moreover, if p = f 1, . . . , is the 
path from p to Vr, then for a we have m{vi) = . . . = m{vr) = 1 while for r we 
have m{vi) = . . . = m{vr) = —1. Writing Ui,Wi for the two siblings of Vi for 
2 < i < k,we see that for all i, for both a and r we have m{ui) 7^ m{vi). Note 
that these events are independent for different values of i. We therefore obtain 
that the probability that Wi, . . . , is such a path is bounded by (2e + 8e^)^~^ 
Since there are 3*" such paths and since P[cr — > r] < 3"'' we obtain that the 
right term of (fTHj) is bounded below by 



Note that the proof above easily extends to the (i-regular tree for d > 3. 
A similar proof also applies to the binary tree T, where g is now defined as 
follows. Look at Tj. for even k. For the boundary vertices define = a^. 
For each vertex v at distance 2 from the boundary, choose three leaves on 



or 




(18) 



(2e + 8e2)i-^>n^('^). 



□ 
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the boundary below it f i, f2, and let be the majority of the values m^.. 
Now continue recursively. 

Repeating the above proof, and letting pi = P[my ^ a^] for a vertex at 
distance 2£, we derive the following recursion: pi < 3(2e)^ + 6(2e)p£_i + 3p^_p 
We then continue in exactly the same way as for the ternary tree. 

4 High temperatures 

Proof of Theorem \1.4\ part Our analysis uses a comparison to block dy- 
namics. 

Block dynamics. We view our tree T = tJ^^ as a part of a larger 6-ary 
tree T^, of height r + 2h, where the root p of T is at level /i in T^,. For each 
vertex f of T*, consider the subtree of height h rooted at v. A block is 
by definition the intersection of T with such a subtree. Each block has a 
rate 1 Poisson clock and whenever the clock rings we erase all the spins of 
vertices belonging to the block, and put new spins in, according to the Gibbs 
distribution conditional on the spins in the rest of T. 

Discrete dynamics: In order to be consistent with we will first analyze 
the corresponding discrete time dynamics: at each step of the block dynam- 
ics, pick a block at random, erase all the spins of vertices belonging to the 
block, and put new spins in, according to the Gibbs distribution conditional 
on the spins in the rest of T. 

A coupling analysis. We use a weighted Hamming metric on configu- 
rations, 

c?(a,r/)=5^AHl(a,^r/„), 

V 

where |f | denotes the distance from vertex v to the root. Let = 1 — 2e 
and A = 1/Vb. Note that bX6 < 1 and 6 < X. Starting from two distinct 
configurations a and rj, our coupling always picks the same block in a and 
in 7] and chooses the coupling between the two block moves which minimizes 

We use path-coupling i.e., we will prove that for every pair of configu- 
rations which differ by a single spin, applying one step of the block dynamics 
will reduce the expected distance between the two configurations. 

Let V be the single vertex, such that ^ rjy. Then d{a, rj) = A'^L Let B 
denote the chosen block, and o"', rj' be the configurations after the move. In 
order to understand {(y',ri'), we will need the following Lemma. 
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Lemma 4.1. Let T be a finite tree and let v ^ w be vertices in T. Let 
{Pe > 0}eeE{T) be the (ferromagnetic) interactions on T, and let {— oo < 
H{u) < oo}u(^v{T) be an external field on the vertices of T . we consider the 
following conditional Gibbs measures: 

^+,h: The Gibbs measure with external field H conditioned on = 1. 
fj^-^H- The Gibbs measure with external field H conditioned on = —1. 
Then, the function fi^^^Ww] — achieves its maximum at H = 0. 

Before proving the Lemma, we utilize it to prove Theorem ll.4t part 121 
There are four situations to consider. 

Case 1. if B contains neither v nor any vertex adjacent to v, then d{a', rj') = 
d{a, T]). 

Case 2. If B contains v, then a' = v( and d{a',r]') = = d{a,r]) — A'^L 
There are h such blocks, corresponding to the h ancestors of f at 1,2, ... ,h 
generations above v. (Note that this holds even when v is the root of T or a 
leaf of T, because of our definition of blocks). 

Case 3. If B is rooted at one of f 's children, then the conditional probabil- 
ities given the outer boundaries of B are not the same since one block has 
+ 1 above it and the other block has —1 above it. However both blocks have 
their leaves adjacent to the same boundary configuration. When considering 
the process on the block, the influence of the boundary configuration can be 
counted as altering the external field. Since a and rj have the same external 
fields and the same boundary configuration on all of the boundary vertices 
except V, by Lemma l4.ll conditioning on this lower boundary can only re- 
duce d{a' ,rj'). Therefore, we bound d{a' ,rj') by studying the case where one 
block is conditioned to having a -fl adjacent to the root, the other block 
is conditioned to having a —1 adjacent to the root, and no external field or 
boundary conditions. Then the block is simply filled in a top-down manner, 
every edge is faithful {i.e. the spin of the current vertex equals the spin of 
its parent) with probability 9 and cuts information (the spin of the current 
vertex is a new random spin) with probability 1 — 6. Coupling these choices 
for corresponding edges for a and for rj, we see that the distance between a' 
and rj' will be equal to the weight of the cluster containing v, in expectation 
'Y^. \\'"\~^^b^9^ < A'^' /(I — b\9). There are b such blocks, corresponding to the 
b children of v. 

Case 4. If B is rooted at f 's ancestor exactly h + 1 generations above v, 
then the conditional probabilities are not the same since one block has a leaf 
V adjacent to a +1 and the other block has a leaf adjacent to a —1. There 
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is exactly one such block. Again we appeal to Lemma (4.11 to show that the 
expected distance is dominated by the size of the 6 cluster of w. The expected 
weight of f 's cluster is bounded by summing over the ancestors w of v: 

w j 

i-h\e 



{i - ex-^){i - hxey 

Overall, the expected change in distance is 
E(ciK,,')-ci(<T,,,))< 



i-hxe {i-ex-^){i-hxe) j n + h-i 

If the block height h is a sufficiently large constant, we get that for some 
positive constant c, 

E{d{a', i) - d{a, r/)) < < —d{a, r^). (19) 

n n 

Note that max (i(o",r7) = J2j<r^^'' — V^- Therefore, by a path-coupling 
argument (see ^j) we obtain a mixing time of at most O(nlogn) for the 
blocks dynamics. 

Spectral gap of discrete time block dynamics. The (1 — c/n) con- 
traction at each step of the coupling implies, by an argument from [5 which 
we now recall, that the spectral gap of the block dynamics is at least c/n. 
Indeed, let A2 be the second largest eigenvalue in absolute value, and / an 
eigenvector for A2. Let M = sup^,^ |/(<^) " f{v)\/d'{o', rf) and denote by P the 
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transition operator. Then 



|P/(a)-P/(r/)| . 
A2 iW = sup r Since / eigenvector tor A2 

a,v d{a, 7]) 

< sup g P [(., ,) ^ (. , , )] ^ 

<sup5^P[(a,77)-(a',V)]M-^^^''''') 
= Msup 



a,r, d{a, 7]) 



<{l-c/n)M by P|). 

Thus IA2IM < (1 — c/n)M, whence the (discrete time) block dynamics has 
relaxation time at most 0{n). 

Relaxation time for continuous time block dynamics. The contin- 
uous time dynamics is n times faster than the discrete time dynamics. This 
is true because the transition matrix for the discrete dynamics is M = 1 + 
where I is the 2"-dimensional unit matrix. Therefore 



r2 (block dynamics) = 0(1). 

Relaxation time for single-site dynamics. Since each block update 
can be simulated by doing a constant number of single-site updates inside the 
block, and each tree vertex only belongs to a bounded number of blocks, it 
follows from proposition 3.4 of [25^ that the relaxation time of the single-site 
Glauber dynamics is also 0(1). □ 

Proof of Lemma \4-l\ Reduction from trees to paths. We first claim 
that it suffices to prove the lemma when the tree T consists of a path v = 
vi, . . . ,Vk = w. (see FigureEl). To see this, let Ti, T2, . . . , be the connected 
components of T when the edges in the path vi,V2, ■ ■ ■ ,Vk are erased, s.t. 
f j G Tj for z = 1, 2, . . . , k. Let a be a configuration on f 1, . . . , Vk, and for a 
subgraph J let S{J) be the space of configurations on J. The probabihty of 
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a configuration a on vi, . . . ,Vk is 

^ pvn I 



^ exp ( ) ■ n I Yl exp (7i(T U (7^.)) 

\i=l / i=l \TeS{l 

^ /k-l k ^ 



for some external field {H'^} depending only on and {/3e}, where Z and 
Z' are partition functions and denotes the Hamiltonian. 



Vl 
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Figure 9: Reduction from trees to paths. 
We will now prove the lemma by induction on the length of the path 

Vl, ■ ■ ..Vk. 

Paths of length 2. Assume k = 2. Writing f3 for the strength of {vi,V2) 
interaction, H for external field at w = f2. Then, 

= tanh{(3 + H)-tanh{H - (3). 
It therefore suffices to prove that for > 0, the function 

H ^ g(P, H) = tanh(// + /3) - tanh(i/ - /3) 
has a unique maximum at = 0. Consider the partial derivative, 

gH{(3,H)^ cosh-2 (H + (3) - cosh'^ {H - 13). (20) 

Therefore, if /3 > and > then guiP, //) < and if /3 > and < 
then gniP, H) > 0. Thus = is the unique maximum and the claim for 
k — 2 follows. 
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Induction step. We assume that the claim is true for — 1 and prove 
it for k. We denote v' = Vk-i, ^ = fini-Wv' = 1] and similarly j^. Now, 

= ^if^+,H - fi-,H)Wv']ifJ''+,H - (21) 

Since by the induction hypothesis both multipliers in ()2ip achieve their max- 
imums at = 0, we get that fJ^+^nl^'w] — fJ'-,H[o'w] also achieves its maximum 
at H = 0. ' ' □ 



5 Proof of Theorem 11.51 

Recall that we denoted by cr.,. the configuration on all vertices at distance 
exactly r from p. Also recall that /i is the Gibbs measure which is stationary 
for the Glauber dynamics. We abbreviate J f dfi as n{f )- 

Mutual information and estimates. For Markov chains such as 
{(Jr}, it is generally known L36j that (0) follows from (0)), which itself, is a 
consequence of the following stronger statement: 

There exists c=k > such that for any vertex set A C Gr/2 and any 
functions f,g with fi{f) = ji{g) = 0, we have 

Kf9)<e-'''iKfM9')y^\ (22) 

provided that /(a) depends only on cr^ and g{a) depends only on ar- 
will follow from a more general proposition below. For a set A of vertices in 
a graph G we write diA for the set of vertices f in A for which there exists 
an edge (f , u) with u ^ A. 

Proposition 5.1. Let G be a finite graph, and let A and B he sets of vertices 
in G. Let d he the distance hetween A and B and let A he the maximum degree 
in G. For < c < 1, let 

/(c) = c-logc- 1. (23) 
Let c* he the unique < c < 1 satisfying I{c) = log A and for < c < c* , let 

C(c,A) = (l-e'°sA-/(c))-V2_ (24) 
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Further, let X2 be the absolute value of the second eigenvalue of the gen- 
erator of the Glauber dynamics on G, i.e. A2 = Let f = f{a) depend 
only on the values of the configuration in A and g = g{a) depend only on the 
values of the configuration in B. If fi{f) = fi{g) = 0, then 

^lifg) < l^e-^''- + 2G{c, A)^\d,A\ e'^(i°gA-/(c))^ WfhWgh (25) 
In particular (by letting c = e~^°sA-7-2j 7 > 0, 
Kfg) < (e-'^^^^-P^-'^sA-T-s) + 4v/exp(-(7 + l)d) \d,A\^ WfhWgh- (26) 

Proof of Theorem ] 1.^ Note that \diA\ < \A\ < A*"/^. Therefore, to prove 
dUl) we use (j2ni) with B={v : d{v, 0) =r},d = r/2 and 7 s.t. e"^ > A. □ 

Proof of Proposition 15. 11 We use a couphng argument. Let /x be the Gibbs 
measure on G, and let Xq be chosen according to /i. Let Xt and Yt be defined 
as follows: Set Yq = Xq. For t > 0, let Xt and Yf evolve according to the 
dynamics with the following graphical representation: Each v E G has a 
Poisson clock. Assume the clock at v rang at time t, and let Xt_ and Yt_ be 
the configurations just before time t. At time t we do the following: 

1. If 1; G -B then X„ updates according to the Gibbs measure, and Y^ does 
not change. 

2. \{ V ^ B and Xt-{w) = Yt-{w) for every neighbor w of v, then both X 
and Y update according to the Gibbs measure so that Xt{v) = Yt{v). 

3. li V ^ B and there exists a neighbor w of f s.t. Xt-{w) 7^ Yt-{w) then 
both X and Y update according to the Gibbs measure, but this time 
independently of each other. 

For a vertex v E B we define ty to be the first time the Poisson clock at 
V rang. For any v E G\B, we define ty to be the first time the Poisson clock 
at V rang after min(^^^)gEg t^. Note that Xt{v) = Yt{v) at any time t < ty, 
and that ty depends only on the Poisson clocks, and is independent of the 
initial configuration Xq. We let tA = ^ii^veAtv 

Let J-'t denote the (a-algebra of the) Poisson clocks at the vertices up to 
timet Let (P*/)(a) = E{f{Xt)\Xo = a,J't] and let (g*/)(a) = E[fiYt)\Xo = 
a,Tt]. Also, let (P*/)(a) = E[f{Xt)\Xo = a] and {Q'f){a) = E[/(r,)|Xo = 
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a]. Since for all t the process Yt is at the stationary distribution and 1^1^ = 
Xq\b for all t, we get 

fi[gf] = E[g{Yt)f{Yt)] = E[g{X,)f{Y')] = E[gQ'f]. (27) 

If t < t^, then clearly Xt = Ft on A Therefore, \\{Q^f-P^f)-lt<tAl = 0- 
On the other hand, ||Q*/||2 < \\fh - a.s. and ||P*/||2 < \\fh - a.s. 
This is because the operators / Q^{f) and / — >• -P*(/) given J^t are Markov 
operators and hence contractions. Therefore 

P'f{cr)]'\t < tA) 
P'f{a)]'\t>tA,Xo = a) 

where the first inequality is because Q*/(Xo) — fi^Xo) is a conditional 
expectation of Q^fi^X^) — P^f^Xo), and the second inequahty is because 
{t > Ia} is jFi-measurable. Therefore, by the Cauchy-Schwartz inequality, 

E[(gV - P'm < 2v/P[tA<t]||/||2 Wgh- (28) 

Since 

E[^?PV]<e-^^*||/||2||^||2, 
We infer that from (j2HI) and (j27|) that 

ftlfg] < (e-'^* + 2^P[tA < t]) ||/||2|kl|2 . (29) 

It remains to bound the two terms in the right-hand side of ()29|). 

For < c < c*, we take t = cd. We obtain that the first term is e"'^'^'^^^ 
as desired. It remains to bound P[tA < t]. We note that < ^ only if 
there is some self-avoiding path (sometimes referred to as "path of disagree- 
ment" ) between the A and B along which the discrepancy between the two 
distributions has been conveyed in time less than t. 

Time-reversing the process, this means that first-passage-percolation with 
rate-1 exponential passage times starting at A needs to arrive at distance d 



WQ^f-p^jWl = E([gV™-PV™f) 

< E ([QV(Xo) - P'fiXo)]') 

= Pit < tA) j dii{a)E {[Q'f{a) - 

+ P(t > tA) j rf/i(or)E {[Q'f{a) - 

< 4P[t^<t]||/||^ 
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within time cd. There are at most \diA\ paths of length k for the first- 
passage-percolation for each k > d. Let r(t>,w) be the time needed to cross 
the edge (f , w). For each path fi, t>2, . . . , ffe, 

P {t{vi,V2) + t{v2, v^) + ... T{vk-i,Vk) < cd) < e'''^^"'^ 

where J(c) = c — logc — 1 is the large deviation rate function for the expo- 
nential distribution. Therefore, 

oo 

P{tA<t) < |9,A|5^exp[fc(-/(c') + logA)] 

k=d 

□ 

Plugging this bound into (j29j) . we obtain (j25j) as needed. 



6 Open Problems 

In this section we specify some relevant problems that are still open. 

Problem 1. What is the relaxation time T2{n, b, of the Glauber dynam- 

ics of the Ising model on the b-ary tree of depth n at the critical temperature 

Using the sum of spins as a test function, we learn that Q{logn) is a lower 
bound for r2(n, 6, We conjecture that the relaxation time is of order 

0(log?T,). A weaker conjecture is that 

\og{T2{n,b,b-'^')) ^ ^ 

n— >oo fl 

Problem 2. Fix b, and let 

log(r2(n,6,/3)) 



r2(/3) = lim 



n 



Theorem \1.J\ part [7] tells us that T2{(3) exists and is finite for all [3. Show 
that T2{P) is a monotone function of p. This question is a special case of a 
more general monotonicity conjecture due to the fourth author, described in 
VdU^ . See VJUl where a monotonicity result is proven for the Ising model on 
the cycle. 
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Problem 3. For the Ising model (with free boundary conditions and no exter- 
nal field) on a general graph of bounded degree, does the converse of Theorem 
M.^ hold. i.e., does uniform exponential decay of point-to-set correlations im- 
ply a uniform spectral gap? 

(As pointed out by F. Martinelli (personal communication), the converse fails 
in certain lattices if plus boundary conditions are allowed). 

Problem 4. Recall the general upper bound ne^^^^'^^^'^^^'^ on the relaxation 
time of Glauber dynamics in terms of cut-width from proposition li.il For 
which graphs does a similar lower bound of the form T2 > e^^'^^^ (for some 
constant c > 0) hold at low temperature ? 

Such a lower bound is known to hold for boxes in a Euclidean lattice, 
our results imply its validity for regular trees, and we can also verify it for 
expander graphs. A specific class of graphs which could be considered here 
are the metric balls around a specific vertex in an infinite graph T that has 
critical probability pd^) < 1 for bond percolation. 

Remark. After the results presented here were described in the extended 
abstract [19j, striking further results on this topic were obtained by F. Mar- 
tinelli, A. Sinclair, and D. Weitz |2n|- For the Ising model on regular trees, in 
the temperatures where we show the Glauber dynamics has a uniform spec- 
tral gap, they show it satisfies a uniform log-Sobolev inequality; moreover, 
they study in depth the effects of external fields and boundary conditions. 
Acknowledgment. We are grateful to David Aldous, David Levin, Laurent 
SalofF-Coste and Peter Winkler for useful discussions. We thank Dror Weitz 
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